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ABSTRACT 



Isospin violation in the p(770) mass and width is considered within the S matrix approach using 
combined fits to the e + e~ tt + tt~ and t~ — > Vt1t~tt data performed by the ALEPH collaboration. 
We show that the pole position following from the parameters obtained from the ALEPH fits are 
not sensitive to the details of the parametrization. In this context, we have found that the pole 
mass difference and the pole width difference between the charged and neutral p are consistent with 
zero. We show that a one loop calculation including vector, axial vector and pseudo-scalar mesons 
can satisfactorily describe the observed isospin breaking. We also give an estimate for the mass 
difference between the neutral and charged states of the ai(1260). 
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I. INTRODUCTION 



Isospin is in very good approximation a symmetry of the strong interactions. At the fundamental level, isospin 
breaking occurs as a consequence of the electromagnetic interactions and the mass difference of the u and d quarks. 
The numerical value of the fine structure constant and the fact that \m u — rrid\ is small compared to Aqcd explains 
why this is a symmetry good at the few percent level. In this work we are concerned with the isospin and its breaking 
in the p resonances. The experimental tests of this symmetry requires the precise determination of the mass and width 
of the charged and neutral states, a non trivial task since a precise definition of these quantities must be adopted. 

The values compilated by the Particle Data Group for the p(770) mass and width are spread over a wide range. 
This is an unpleasent characteristic since one expects these intrinsic properties to be independent of the process where 
the resonance is observed and also to be model and background independent. This is related to the fact that different 
expressions are used to fit the data leading to different values for the parameters - that one is tempted to identify 
with the physical mass and width - which however are relevant only for the model under consideration. Although this 
is a practical way to deal with the p meson physics, i.e. use a given parametrization and the corresponding values for 
the mass and width, the process and model independence of the physical quantities associated to any resonance are 
relevant from the conceptual point of view. 

The analysis of the experimental data is based on theoretical assumptions such as a parametrization derived 
in terms of analyticity and unitarity Q , a Breit-Wigner which incorporates an energy dependent width reflecting the 
p-wave nature of the p and higher isovector- vector resonances || or still phenomenological effective field theories [Q . 
There is a priori no reason to favour any of these approaches in the description of physical process, on the contrary one 
would expect them to yield complementary or at least consistent information. All of them share two characteristics: 
i) they involve a momentum dependent width T(s), and ii) they define the resonance mass M as the energy for 
which the real part of the propagator is zero and the width as T(s = M 2 ). We will refer collectively to them as the 
conventional, as opposed to the pole approach to be discussed next. On the other hand, the S matrix formalism Q 
relates the mass and width of the resonance to pole(s) of the amplitude in the complex plane. The pole associated to 
a given resonance is independent of the process where it is observed [||^] , thus this scheme offers the possibility to get 
an unambiguous determination of the intrinsic properties of the resonance. For the p(770) this has been done in Jl0| . 

The S matrix formalism not only provides a convenient parametrization in terms of which the electromagnetic form 
factor of the pion can be expressed, but it can also be applied to any of results of the conventional approach. In fact 
we have two possibilities: the first is to make an expansion around the pole position of the expressions used in the 
conventional approach and then fit the data. We already know the result of such a procedure since in all cases the 
form factor will reduce to the S matrix formalism, and the fit using this scheme was performed in JTcj] . Alternatively, 
fits to the data can be carried using the conventional parametrization and then, given the expression used and the 
parameters obtained in fitting the data, one can determine the associated pole. It is our purpose in this letter to 
perform the latter analysis considering fits reported for e + e~ — ► ir + ir~ and for the tau decay r~ — ► v T it~i!p which 
allow us to make a statement about the size of isospin violation in the p^- p° system within a framework in which 
the masses and widths are unambiguously defined. 

We also consider the problem from the theoretical point of view. To this end we take into account both possible 
sources of isospin violation: the m u — mass difference and the electromagnetic interactions, which is implemented 
at the hadron level by the use of vector meson dominance. The former contribution is known to vanish to lowest order 
implying thus a finite result for the latter. The explicit calculation of the one loop electromagnetic self energy shows 
that a finite result can be achieved only when the anomalous magnetic moment k entering the three p vector meson 
vertex (Vector Meson Dominated p — p — 7 vertex) takes the value K = 1. We briefly comment on this point. 



II. M AND T IN THE CONVENTIONAL AND POLE APPROACHES. 



Within the conventional approach, near the p resonance, the electromagnetic form factor of the pion is 
parametrized in terms of a Breit-Wigner with a momentum dependent width: 

F ^ = 7^ ; D ^ = s ~ ^ + l ^p( s ) f p( s ) (!) 
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Although some expressions ocurring in the literature seem to differ from Eq.(|j]) they can be cast in that form by an 
appropiated redefinition of A(s) and T(s). For example if instead of A(s), we use A(s) = A(s)(l — ix(s)), D(s) should 
be changed to: 

D(s) = s - M 2 (s) + M p (s)T p (s)x(s) + i(M p (s)f p (s) - (s - M 2 (s))x(s)). (2) 



Since in the conventional approach the mass and width of the resonance are defined respectively by 

Re{D(m 2 p j) = ; and Im(D(m 2 p ) = m p T p , (3) 

the freedom in the choice of parametrization (compare Eqs.([l]]^)) implies an arbitrariness in the mass and width 
deffinition. On the other hand, the S matrix approach ||||] avoids this problem by considering a simple observable 
pole with s independent width and residue: 

F r (a) = + B(a), (4) 

S Sp 

where s p is the pole position, R p the pole residue and B(s) is the remaining part around the pole. In order that this 
description makes sense the remaining part -B(s) (which is fixed by the fit to the data) should be a soft function of s 
affecting minimally the pole position, except for obvious resonance effects as in p — uj mixing JT(| . In this formalism 
the mass M and width T of the resonance are deffined through: 

Sp = M 2 - iMT. (5) 

The fact that the background is a soft function arround the pole (in this case for ^fs ?^ip), does not imply that 
it can not have poles at other energies, however these can be included in an appropiated way in the background. 
Different parametrizations within the conventional approach correspond to different backgrounds in the S matrix 
formalism. Thus even though the mass and width values obtained from a given "conventional" treatment are only 
relevant for that particular parametrization, on physical grounds one expects that all of them have a pole in the same 
position providing thus, through Eq.(|^), the physical mass and width of the resonance. Notice that in finding the 
pole D(s p ) = no information about the overall normalization nor about A(s) is required, as one would expect from 
the <S* matrix approach. 



The cleanest determination of the p mass and width comes from the e + e anihilation and tau-lepton decay data, 
which are in agreement with each other up to an overall normalization factor @. Although several analysis for the 
e + e~ — * 7r + 7r~ ||,||,0,[ll] lj] and t~ — > is T ir~ir° [^,|| data have been performed, we concentrate in the results of the 
ALEPH collaboration |2] which seems to us the best suited for a study of isospin breaking since a combined fit to 
both sets of data was carried by these authors. Notice however that the agreement is excellent with the fit to the data 
of the high statistics experiment by the CLEO collaboration [|| . The p^ , p° masses and widths resulting from the 
analysis in [|| and the values associated to the corresponding pole position are summarized in Table 1 for the masses 
and Table 2 for the widths. 



Fit 


m p ±(MeV) 


m p o (MeV) 


m p ± — m p o (MeV) 


conven. 


pole 


conven. 


pole 


conven 


pole 


KS(A = 1) 
GS(A = 1) 
GS(A = 0.45 ±0.11) 


773.4 ±0.9 

775.7 ±0.9 

783.8 ±3.0 


757.0 ± 1.3 

758.1 ±1.3 
758.3 ±5.4 


773.4 ±0.7 

775.7 ±0.7 

783.8 ±3.0 


756.9 ± 1.0 

758.0 ± 1.0 

758.1 ±5.4 


0.0 ± 1.0 
0.0 ± 1.0 
0.0 ± 1.2 


0.1 ±1.6 
0.1 ±1.6 
0.2 ±7.6 



Table 1. Masses of the p and p° resonances as obtained by Barate et al J2| (second and fourth column) from a 
combined fit to e + e~ — > -k + tt~ and r _ — » ^ T 7r _ 7r° data. For comparison in the third and fifth column we quote the 
values associated to the corresponding pole. Details about the notation can be found in the quoted literature. 
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Fit 


I>(MeV) 


T p0 (MeV) 


T p ±-T p o(MeV) 


conven. 


pole 


conven. 


pole 


conven 


pole 


KS(A = 1) 
GS(A = 1) 
GS(A = 0.45 ±0.11) 


147.7 ± 1.6 

150.8 ± 1.7 
162.0 ±5.3 


143.2 ± 1.5 

145.3 ± 1.5 
145.1 ±6.3 


147.3 ± 1.3 
150.8 ± 1.3 

162.4 ±5.0 


142.7 ± 1.2 
145.2 ± 1.2 
145.2 ±6.5 


0.4 ± 1.0 
0.0 ±2.0 
-0.4 ±2.5 


0.5 ± 1.9 
0.1 ±1.9 
-0.1 ±8.7 



Table 2. Widths of the p ± and p° resonances as obtained by Barate et al || (second and fourth column) from a 
combined fit to e + e~ — > tv + tv~ and r~ — > i/ t tt~tv° data. For comparison in the third and fifth column we quote the 
values associated to the corresponding pole. Details about the notation can be found in the quoted literature. 



The results are remarkable. While in the conventional approach the masses are spread over a 10 MeV range, the S 
matrix formalism leads values that differ at most in 1.3 MeV. Similarly for the width where the two approaches yield 
values differing, within each scheme, in 14 and 2 MeVs respectively. We can safely conclude that in contradistinction 
to the conventional approach the pole position is independent of the details of the parametrization used to fit the data 
providing thus an unambiguous definition for the mass and width of the p. It is worth remarking that the mass and 
width values we obtained applying the S matrix formalism to different parametrizations of the conventional approach 
agree with previous anlysis performed fully in the S matrix formalism JToj , |]~3| — |l5|| but are systematically below the 
values obtained in the conventional approach and the value reported by the PDG jl| . 

As far as the isospin breaking is concerned, both schemes lead similar results. Although in the conventional approach 
there is a strong dependence upon the parametrization used to fit the data, it is the same for the charged and neutral 
channels and it cancels out when differences are considered. Let us emphasize again that the advantage of the pole 
approach is that it is based on an unambiguous definition of the mass and width and so, besides a nearly vanishing 
p — p° mass difference we have a precise knowledge of each mass and width. The isospin violation in the coupling 
constants, g p , implied by the pole decay rates is also consistent with zero, g 2 p ± and g p0 are equal up to one percent. 

When applying the S matrix formalism to one of the parametrizations used in the conventional approach, as we have 
done, the uncertainty in the pole position follow from the result of the original fit. In order to calculate the errors 
quoted in Tables 1 and 2, we used the following procedure: we obtain the pole position using the central value of 
the parameters obtained from the fit, we then determine the variation of the pole position when the parameters of 
the fit are changed by one standard deviation. We quote for the pole position error the maximun variation found 
in this way. One further comment regards the A parameter [ fl2| used in the second version of the Gounaris-Sakurai 
parametrization A modify the s dependence of the p decay width and it also changes the analytical structure 
obtained in an effective field theory. Notice that the errors are much larger when A is left as a free parameter, i.e. 
when the s dependence is not the one predicted by the effective field theory . 



III. M p ± - M p0 MASS DIFFERENCE WITHIN A VECTOR MESON DOMINANCE MODEL. 

In this section we consider the theoretical evaluation of the p ± — p° mass difference. In order to calculate the 
contribution due to the m u — m d mass difference, we need to evaluate the matrix element of the flavour symmetry 
breaking Hamiltonian Hb = m u uu + m d dd + m s ss. The matrix elements between p states can be written as: 

,/ _ T, _ m 1 / s .uu — dd. 
< p\(m u uu + m d dd + m s ss)\p > = -j={rn u - m d ) < p\ -j= \p > 

1 . . .uu + dd— 2ss . 1 . . .uu + dd + ss. , . 

+ -j=(m u + m d -2m s ) < p\ -j= \p > +-j=(m u + m d + m s ) < p\ -j= \p > . (6) 

The operators entering in the two last matrix elements are singlet under SU(2), therefore they give the same con- 
tribution to the p and p° and thus, do not contribute to the mass difference. On the other hand the first matrix 
element vanishes identical as consequence of G parity. We therefore conclude that to lowest order in (m u — m d ) the 
p ± — p° mas difference vanishes: 
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m p ± - m p o\ mu jk md = 0. 



(7) 



This result has an important implication. Since the m u ,d quark masses have completely disappeared from our 
calculations, there are no fundamental constants of the underlying QCD theory which could be used to renormalize, 
i.e. to absorbe possible divergencies appearing in the calculation. Therefore the electromagnetic self-mass contribution 
must be finite |L6| |. Later on we will see that the same conclusion can and has been reached ]T7[ using different 
arguments. 

Let us now consider the second source of isospin breaking, the electromagnetic interactions. The lowest order 
contributions are shown in Fig. l(a— e), where P stands for a pseudoscalar meson which in our case can be either it, 
rj or 77/, ai denotes an axial vector and V = p, ui or <j> enters through Vector Meson Dominance. The first diagramm 
is non vanishing for the charged p and vanishes for the p . Diagram b is non vanishing (V ~ lu or <fi) both for the 
charged and neutral contribution. Notice furthermore that the r\ or 77/ can be intermediate state only for the neutral 
p° self-energy. We will neglect the contribution of the 77 and 77/ as we expect these to be small (one hundred times 
smaller than the contribution of diagram e). On the other hand, the contribution from the loop with an internal pion 
yields the same contribution for the charged and neutral p's. The third diagramm vanishes identically when V = p 
due to G parity and gives the same contribution to the p and p° self-mass when V = u> or <f> and therefore does 
not contribute to the mass difference. Figures d and e are non-vanishing and contribute only to the charged and 
neutral meson mass respectively. Thus, in the model we are considering, the p - p° mass difference is given by the 
contributions shown in Fig. 1 (a, d and e). To describe the photon p vertex we use Vector Meson Dominance, so that 
the p form factor is modeled by a p propagator. The Feynman rules for the electromagnetic interaction of a charged 
massive spin S = 1 vector meson can be found in [ ^8[ . These rules depend upon the anomalous magnetic moment 
k of the vector meson. Using these rules, we find that the correction to the mass is given by (for later use we 
introduce obvious notation for the loop integrals of Ti,T% and T^k a kp): 





P P 

JXTUIAAAAT 

e) 
Fig. 1 

One loop electromagnetic contribution to the self-mass of the p. Vector Meson Dominance is used to describe the coupling of 
the photon to the hadrons. P stands for a pseudoscalar (ir, rj or r)l), V for a vector meson (p, to or cf>) and a for the axial 

vector ai(1260). See main text for further details. 
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id 



7 _ 



e 2 M 4 e°'(p)e* p (p) J -^-^ (Ti g a p + T 2 g a p + T 3 k a kp) 
aM 4 



d 4 k 



An 



-(h + t 2 + t 3 ) 



(8) 



where: 



Ti = 



Am 2 



k 2 [k 2 — M 2 ] 2 [k 2 + 2pk] k 2 [k 2 - M 2 ] 2 [k 2 — M 2 ] 2 [k 2 + 2pk] ' 



(9) 



T 2 = 



K 2 k 2 



3k 2 - 8k + 4 

[k 2 - M 2 } 2 [k 2 + 2pk] Am 2 [k 2 - M 2 } 2 [k 2 + 2pk\ Am 2 [k 2 - M 2 } 2 



n = 



k 2 + An - A 



k 2 - 2k + 1 



k 2 [k 2 - M 2 ] 2 [k 2 + 2pk] m 2 [k 2 - M 2 ] 2 [k 2 + 2pk] ' 



Remark that the T± contribution to the integral is exactly the same entering in the calculation of the electromagnetic 
mass difference of the pion Jig ). Notice furthermore that as it stands, the expression involves ultraviolet divergences. 
For complctness we quote the analytical expression for (5 m 2 as a function of n and the p (m) and p°(M) masses: 



2 

M 2 



1 , M 2 
2m z m z 



M 2 +2m 2 _ x R 

,R tan ( j^ ), 



! M 4 



(10) 



A k 2 

t 2 = t(k 2 - 2k + 1) - — j 

m Am* 



M 2 M 2 

2+3 T )ln^r 

m z m- 2 



2(M 2 ~ Am 2 ){M 2 - m 2 ) , R 



2 R 



t 3 



(k 2 -2k+ 1) 

4m 2 
(k 2 +Ak- A) 
Am 2 



3 + A 



M 4 - Am 2 M 2 



2m 4 

M 2 -3m 2 M 2 
3m 2 m 2 



2m 4 , M 2 
In — o 



M 2 

-.2 



m 

2 J\/f2V/l^ 2 



2m 2 - M 2 „ _-, . R . 
M 4 R Un ~ W 



2{m 2 - M 2 )(Am 2 - M 2 ) R . 

-ZNPR tm { W 2) 



with: 



2 M 2 

7 + ln(An) — In — 5- 

— n [i 2 



R= V /Af 2 (4m 2 - M 2 ). 



(11) 



The self-energy involves ultraviolet divergencies - contained in the A factors - which are absent for k — 1. On 
the other hand the finiteness of the p electromagnetic self-mass is required not only from the reasoning presented 
paragraphs above. In Ref. |l7j the same conclusion was reached using completely different arguments (based on 
dispersion relations and Regge pole theory). Furthermore one expects the p — p — p vertex to be symmetric in the 
indices, which is achieved only if k = 1. Taking that value for k and m = M = 757.5 MeV we get: 



51: 



2m,r 



am p 

~8^r 



2 + irV3- 



1.49 MeV. 



(12) 



G 



On the other hand, in the model we are considering the correction to the p° meson mass is given by the diagramm 
in Fig. 1(e). which leads to: 



2 it a m, 



^. - 2^ - --^rr (13) 

With 4.83 < f p < 4.91 (as determined in Ref.( |^3))) from fits to e + e~ — * tt + it~ data in the p region in the 
framework of the S matrix formalism, we use the central value f p — 4.87 so that putting all toghetcr we finally obtain: 

Am p = m p ± - m p o = S m ± - S m = 0.02 ± 0.02. (14) 



Our result is consistent with a previous calculations by Bijnens and Gosdzinsky [jl9| who using a different approach 
find -0.7MeV < m p ± - m p o < OAMeV. 

A similar analysis can be carried for the a± (1260). To obtain an estimate we neglect the analog of diagram in Fig. 
lb which contributes to the self mass of the charged but not of the neutral ai(1260). G parity requires V = p and 
the direct coupling of the photon to the neutral a\ (Fig. le) is not alllowed. The contribution to the ai(1260) mass 
difference from diagrams la and Id is obtained replacing m p by m a in Eq. (|l2|), which leads to: 

Am Q = Sm a ± - Sm a o = 5m a ± = 2A2MeV. (15) 



IV. SUMMARY 



In this work we considered the difference in masses and widths of the charged and neutral p(770) mesons. The S 
matrix formalism provides a framework where these intrinsic properties are unambiguously defined. The determination 
of the mass and width of the p in that framework have been presented in Ref. [ p~0|Jl3| . Here we go one step further 
by taking as starting point the combined fit to the e + e~ — ► tt + tt~ and r~ — > ^ r 7r - 7r r data performed by the ALEPH 
collaboration 0]. The anlaysis in |2| was carried using the Khun-Santamaria formalism and two versions of the 
Gounaris-Sakurai H parametrization of the pion form factor. We have shown that the pole position following from 
the parameters obtained from the fit to the data are not sensitive to the details of the parametrization. In this 
context, we have found that the mass difference between the charged and neutral p is consistent with zero very much 
as the width is. Our results for the p mass and width agree with the pole position found when the pion form factor is 
parametrized from the very begining according to the S matrix formalism and are smaller than the values reported 
by the PDG. 

From the theoretical point of view, we considered the isospin violation induced by the m u — mass difference 
which is known to vanish to the lowest order. We have argued that the electromagnetic self-mass should be finite and 
we have shown that this result is obtained when the elctromagnetic coupling of the p is modeled using vector meson 
dominance and the p anomalous magnetic moment k takes the value k = 1. The model satisfactorily describes the 
observed isospin breaking. 

Let us finally comment that the parametrization independence of the pole position holds when a larger set of fits is 
considered. Indeed when the typical parametrizations used in Ref. p|-^,[7PT2]| are analyzed along the lines of this work 
(except for the fit called VMD2(NU) |7]], which as pointed out by the author leads to a x 2 = 148/77 compared to the 
X 2 ~ 1 of all the others fits) the mass values associated to the pole are spread over a five MeV region arround 758 
MeV, in contrast to the 25 MeV arround 767 MeV obtained in the conventional approach. A similar situation holds 
for the width where the spreads are respectively 7 and 22 MeV and the central values 143 and 151 MeV respectively. 
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